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Abstract: The multi cellular tumor spheroid (MCTS) model has been used for decades with proven superiority over mon- 
olayer cell culture models at recapitulating in vivo tumor growth. Yet its use in high-throughput drug discovery has been 
limited, particularly with image based screening, due to practical and technical hurdles. Here we report a significant ad- 
vance in utilizing live MCTS models for high-content image based drug discovery. Using a validated GFP reporter 
(CK5Pro-GFP) of luminal breast cancer stem cells (CSC), we developed an algorithm to quantify changes in CKSPro- 
GFP expression levels for individual Z-stack planes (local) or as maximal projections of the summed Z-stacks (global) of 
MCTS. From these image sets, we can quantify the cross-sectional area of GFP positive cells, the fluorescence intensity of 
the GFP positive cells, and the percent of spheroid cross-sectional area that expresses CKSPro-GFP.We demonstrate that 
acquiring data in this manner can be done in real time and is statistically robust (Z'=0.85) for use in primary high-content 
screening cancer drug discovery. 

Keywords: Multicellular Tumor Spheroid (MCTS), High-Throughput Screening (HTS), High-Content Screening (HCS), Can- 
cer Stem Cells (CSC), Cytokeratin 5, "Systems Biology", Luminal Breast Cancer, and Cancer Drug Discovery. 



INTRODUCTION 

It has been forty-three years since Sutherland and co- 
workers first described the in vitro 3D tumor model, now 
referred to as the multicellular tumor spheroid (MCTS) [1, 
2]. From this very first report and subsequent articles the 
MCTS model has displayed properties that significantly dif- 
fer from monolayer (2D) cancer cell culture; notably a mor- 
phology that closely resembles in vivo human tumor growth 
and resistance to conventional chemotherapy and radiation 
[1, 3-10]. Accordingly, it is well accepted that the MCTS 
model may be useful in bridging the gap of information that 
exists between 2D-cell culture models and in vivo models 
[11, 12]. Indeed, MCTS models are highly amenable to het- 
erotypic cell culture, and have played an important role in 
characterizing the function of the microenvironment in tu- 
mor progression as well as promoting invasion and metasta- 
sis [13-16]. When compared to 2D-cell culture models, mul- 
tiple studies have demonstrated that 3D cell culture models 
better recapitulate in vivo tumor gene expression profiles [17, 
18]. Furthermore, the MCTS appears to be ideal for studying 
subsets of cells that exist in heterogeneous tumors in vitro, 
particularly those cells that can self renew, initiate tumor 
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growth, resist cell death, and can promote tumor recurrence 
and metastasis, termed herein as cancer stem cells (CSCs) 
[19-21]. Hence, the MCTS may also be an ideal model for 
drug discovery and development targeting CSCs. 

Advances in technology and methods to culture MCTS 
suitable for high-throughput screening (HTS) and more re- 
cently image-based high-content analysis have generated a 
surge of interest in utilizing the MCTS model for novel drug 
discovery. Recently, we reviewed this topic and refer the 
reader to this article for a comprehensive overview [5]. In 
addition to its superiority over 2D-ceIl based models in mim- 
icking in vivo tumor growth, the MCTS model offers a "sys- 
tems biology" [22, 23] approach to cancer drug discovery in 
the context of the tumor microenvironment packaged in one 
assay. For example, the tumor microenvirormient has a role 
in regulating and promoting cell signaling, gene expression, 
and overall tumor morphology including invasion and metas- 
tasis. If we consider the spheroid as a tumor and the extracel- 
lular matrix(ECM) as part of the microenvironment then 
together this could be considered a basic system. Incorporat- 
ing heterotypic cell culture and other components of the mi- 
croenvirormient would expand such systems in the future. 
Heterotypic MCTS culture methods are adaptable to high- 
throughput technologies and screening platforms including 
image-based high-content screening (HCS) [24-28]. Never- 
theless, significant hurdles exist with HCS data acquisition 
and automated analysis, limiting MCTS systems to endpoint 
analysis (e.g. MCTS viability) with no spatiotemporal in- 
formation. Therefore, developing methods to analyze live 
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Fig. (1). Live MCTS expressing GFP reporter of luminal breast CSCs imaged in 3D. Maximal Projections of a CK5Pro-GFP-T47D 
spheroid treated with 100 nM progesterone (P4) over 72 h, resulting in global expansion of the CSC sub-population (top panel). RU486 (RU) 
1 l^lM blocks global expansion of this CSC phenotype in CK5Pro-GFP-T47D spheroids (bottom panel). See supporting information for time 
lapse videos of P4 induced CSC expansion and inhibition using 1 |.lM isotretinoin. 



MCTS models locally (within cells and zones of spheroids) 
and globally (whole spheroids) using automated high-content 
imaging will effectively advance the MCTS model to be 
used as an HCS system for cancer drug discovery. To this 
end, we report a novel MCTS model designed to identify 
small molecule modulators of a sub population of CSCs ex- 
isting within a population of luminal breast cancer cells. This 
model system utilizes cytokeratin 5 (CK5) as a biomarker 
readout and surrogate reporter for this CSC phenotype in 
T47D luminal breast cancer cells. 

CK5 marks bi-lineage progenitor cells in the normal 
breast capable of differentiating into luminal or ba- 
sal/myoepithelial cells [29, 30]. Increased expression of CK5 
is correlated with poor prognosis across all breast cancer 
subtypes [31]. Poor prognostic basal-like breast cancers, the 
major subtype of triple negative breast cancer, have the 
highest CK5^ cell content [32, 33]. Many luminal estrogen 
receptor (ER) ^breast tumors also contain subpopulations of 
CK5^ cells that are enriched in the expression of stem and 
mesenchymal genes, and are more quiescent, self-renewing, 
and drug resistant compared to intratumoral CK5~ cells [34, 
35]. Furthermore, the CK5^ cell population in luminal breast 
cancer expands upon treatment with progesterone (P4), a 
pro-tumorigenic hormone under some contexts [34-37]. Re- 
cently, we validated CK5 expression as a surrogate reporter 
of this CSC population in T47D luminal breast cancer cells 
utilizing the CK5 promoter (CKSPro) enhanced green fluo- 
rescent protein (GFP) reporter stably transduced into T47D 
cells yielding the CK5Pro-GFP-T47D cell line [38]. To vali- 
date this reporter for HCS we utilized DMSO as a negative 
control and RU486, a potent inhibitor of progesterone recep- 
tor (PR), as a positive control. RU486 blocks P4 induced 
expansion and transformation of the CSC phenotype with a 
Z'factor > 0.50 per plate tested. Pilot screening with this 



reporter via a 2D-cell based HCS yielded four hits including 
three retinoids. Retinoic acid proved to be a potent inhibitor 
of CK5 expression (IC50 = 800 pM) and CSC phenotype ex- 
pansion. We have now adapted this reporter for use with 
T47D MCTS in 96-well plates, cultured uniformly, and in- 
corporating an extracellular matrix (ECM). We have devel- 
oped methods to measure and quantitate the modulation of 
this CSC population using retinoids in live MCTS both lo- 
cally and globally, and have validated these methods for 
HCS-based drug discovery. 

MATERIALS AND METHODOLOGY 
Materials 

The CK5Pro-GFP-T47D cell line was produced and 
maintained as previously described [38]. Cell culture rea- 
gents and drug compounds were obtained and used as de- 
scribed below or as directed from the manufacturer, includ- 
ing: Dulbecco's Modified Eagle Medium (DMEM) from 
Coming Life Sciences (Union City, CA, catalogue#15-013- 
CM), phenol free-DMEM from Sigma Aldrich (St. Louis, 
MO, catalogue#D1145), HyClone fetal bovine serum (FBS) 
from Thermo Fisher Scientific (Rochester, NY, cata- 
logue#SH30071), Matrigel™ (Bedford, MA) Progesterone 
(P4) from Sigma Aldrich (St. Louis, MO, catalogue* P8783), 
13-cis-Retinoic acid (isotretinoin) from Sigma Aldrich (St 
Louis, MO, catalogue* R3255), RU486 from Sigma Aldrich 
(St Louis, MO, catalogue* M8046), dimethyl sulfoxide 
(DMSO) and ethanol (EtOH) from Thermo Fisher Scientific 
(Rochester, NY, catalogue* BP231), bisBenzimide Hoechst 
33342 trihydrochloride from Sigma Aldrich (St. Louis, MO 
catalogue* B2261). Unless otherwise stated below, black 96- 
well microplates were used for imaging and were obtained 
from Greiner Bio-One (Monroe, NC, catalogue* 655090). 
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Fig. (2). High content imaging to quantitate tlie lieterogeneity of reporter expression in individual splieroids using all three dimen- 
sions. A) Maximal projection images of spheroids where Z-stack images were taken at steps of 2.6|a.m over a Z span of 240|a.m. Individual 
planes were combined in a maximal projection. High content output parameters shown across each focal plane of spheroids shown in B) The 
sum of fluorescence intensity across each plane, and C) The area of CK5Pro-GFP expression in |a.m". Compound abbreviations are: RU486 
(RU), Isotretinoin (RA), progesterone (P4), and DMSO is the vehicle control. 



CK5Pro-GFP-T47D MCTS Culture 

Single uniform CK5Pro-GFP-T47D breast carcinoma 
spheroids were cultured in 96-well microplates using a modi- 
fied liquid overlay method previously described [6]. Briefly, 
each well of a 96-well flat bottom plate was coated with 
50(a,L of 1.5% agarose solution in serum-free and phenol red- 
free DMEM. During the coating process, the agarose solu- 
tion was maintained at >60 °C followed by cooling and set- 
ting at room temperature for 40 min. CK5Pro-GFP-T47D 
cells were plated at a density of 10,000 cells per well in 
50(xL phenol red-free DMEM supplemented with 10% FBS. 
Next, the plate was centrifuged at 1000 rpm for 15 min at 
room temperature to induce cell aggregation. Aggregated 
cells were coated with 50(a,L of a 10% solution of growth 
factor-reduced MatrigeF"^ (protein concentration: 7.1 



mg/mL) in phenol red-free DMEM supplemented with 10% 
FBS, resulting in a final concentration of 5% MatrigeF"^. 
After two days of incubation, an additional 50|j,L of 5% Ma- 
trigeF"^ containing medium was added to each well to main- 
tain a 5% MatrigeF"^ concentration in a final volume of 
150|j,L per well. Before treatment, spheroids were cultured 
for 7 days yielding an average diameter of 1 mm under stan- 
dard tissue culture conditions (37°C, 5% C02). 

CK5Pro-GFP-T47D MCTS Treatment and Time Lapse- 
Video Capture 

To generate time-lapse videos, spheroids were transferred 
into a glass bottom 96-well plate (Whatman, catalogue# 
7716-2375) and treated with vehicle (EtOH), lOOnM P4 and 
lOOnM P4+l|a,M isotretinoin. Following treatment, spheroid 
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images were captured every 15 min for 12.5 h using a 31 
MARIANAS inverted Spinning Disk microscope. 

CK5Pro-GFP-T47D MCTS Treatment and fflgh-Content 
Imaging 

Spheroids were treated with lOOnM P4 and varying con- 
centrations of RU486, isotretinoin, 5-fluorouraciI (5FU), or 
Tamoxifen, ranging from O.OOlnM to lOjxM, over a 72 h 
timecourse, using DMSO (0.5% final volume [38]) as the 
solvent vehicle. Spheroids were either imaged directly or cell 
nuclei were stained with 10(a,M Hoechst 33342 (bisBen- 
zimide H 33342 trihydrochloride), at 37 °C for 30 min, be- 
fore live spheroid image based analysis. Spheroids were 
transferred to 96 well plates along with 50-100|a,L of media 
and imaged using an Operetta (Perkin Elmer) equipped with 
a live-cell chamber that was pre-equilibrated to 5% CO2 and 
37 °C. For detailed global image acquisition, a lOX High- 
NA objective was used in spinning-disc confocal mode. Z- 
planes were imaged in 2.6^m steps. For whole-plate imag- 
ing, a lOX Long- WD objective was used in spinning-disc 
confocal mode with 8 Z-planes imaged for each field, with 
25|xm steps between Z-planes. 

CK5Pro-GFP MCTS High Content Analysis 

Acquired images were analyzed using Columbus soft- 
ware (Perkin Elmer). Names used for general high content 
output fields are listed in parentheses below. Analysis steps 
include region-based image segmentation based on regions 
of green fluorescence (GFP Region). Selected regions were 
tiled with spots, and intensity and morphological properties 
were calculated for all spots. Spots were selected based on 
spot size (px^) and spot intensity in order to filter out non- 
spheroid related or background selected fluorescence (Se- 
lected GFP Region). High content fields were summed 
across the total area of selected spots to determine the 
"Summed GFP Region: Intensity Sum" and the "Sum GFP 
Region: Area" parameters. For the dose-response and Z- 
prime figures, analysis was done on a maximal projection 
image representing the peak values of each of 8 Z-planes 
projected onto a single image. For the detailed analysis of 
whole spheroids, a similar analysis was done using Harmony 
software (Perkin Elmer), where each Z-plane was analyzed 
individually. 

Statistical Analysis 

Z-primes were calculated using the formula (Zhang 
etal.) [39]: 

Z' Factor = l-'-^^P^ 

Five replicates each of positive controls and negative 
controls were used. Positive controls shown were 10 |iM 
isotretinoin and negative controls included 0.001 nM isot- 
retinoin. For IC50 determinations, each data point is 
representative of data from five spheroids in a single 96-well 
plate. Output values from the analyses above were plotted on 
a semi-log scale and fit to a 4-parameter sigmoidal dose re- 
sponse curve (Graphpad Prism). 



RESULTS AND DISCUSSION 

Although the study of MCTS as a model system for tu- 
mor biology has been pursued for over forty years, more 
recent advances in imaging and screening instrumentation, as 
well as in analytical software, has brought the MCTS to the 
brink of being a valid primary screening tool for cancer drug 
discovery. The advantages to image-based primary screening 
in MCTS are several- fold, and include: the ability to distin- 
guish region specific, local effects of small molecules on 
MCTS phenotype, and whole-spheroid global effects. This 
approach can also provide information on small-molecule 
permeabihty into biologically relevant tissue models as well 
as on the toxicity of screening compoimds at the primary 
screening level. 

The difficulties involved in primary image-based screen- 
ing of MCTS should not be underestimated; however, we 
propose that these are primarily issues of practicality. Here 
we present an initial pathway to large-scale primary screen- 
ing of MCTS with an eye towards multiparametric readouts 
with sub-spheroid spatiotemporal resolution. We first dem- 
onstrate that live MCTS can be imaged in an environmen- 
tally controlled chamber of a high-content imaging system 
(Fig. 1). We chose the CK5Pro-GFP-T47D cell line as the 
MCTS model because of the clinical relevance of CK5 ex- 
pression in breast tumors [38]. 

This reporter of the CSC phenot5^e expresses in only a 
subpopulation of the cells that comprise the individual 
MCTS when treated with P4. We can block P4 induced ex- 
pansion of this CSC phenotype using the potent PR inhibitor 
RU486. In a previous screen, we identified three retinoids 
that do not inhibit PR directly but are also potent inhibitors 
of CK5 expression and expansion of the CSC phenotype 
[38]. Here, we demonstrate the efficacy of one of these reti- 
noids in MCTS, using time lapse videos of P4 and 13-cis- 
Retinoic acid (isotretinoin) treated spheroids over a 12 h 
time course (supplemental material). Individual CKSPro- 
GFP-T47D spheroids were cultured in the presence of 5 % 
ECM (MatrigeF^) to a diameter of ~ 1mm. When P4 is 
added to individual spheroids there is a significant upregula- 
tion of CK5Pro-GFP activity. However, isotretinoin instantly 
blocks P4 induced CK5 expression preventing expansion of 
the CSC phenotype over this 12 h time course. The hetero- 
geneity of CK5 expression in MCTS makes it an excellent 
candidate for image-based analysis, where spatial resolution 
of CK5 expression in live spheroids can be examined and 
quantified in much greater detail than with an endpoint read- 
out such as whole-well luminescence from a luciferase re- 
porter. This will also allow us, when fully validated for 
screening, to separately analyze individual spheroids should 
there be multiple spheroids in a single well of a multiwell 
plate, which will significantly ease the burden of using an 
automated high-throughput spheroid manipulation approach. 

In Fig. (2), we show that our image-based approach is 
capable of not just observing, but also quantifying changes in 
CK5Pro-GFP expression in individual Z-planes using intact, 
live MCTS. The heterogeneity of CK5Pro-GFP expression in 
the X and Y dimensions of each individual plane that we 
observe by eye is recapitulated in the Z-dimension, as shown 
quantitatively in the histograms. Furthermore, we show that 
overall CK5Pro-GFP expression level can be quantified for 
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Fig. (3). Validation testing witli clinically used cancer drugs. To show that CK5Pro-GFP is a selective reporter for the CSC phenotype we 
utihzed the chnically used cancer drugs 5FU and Tamoxifen. MCTS were treated with 100 uM P4 and 1 |lM 5FU or Tamoxifen for 72 h. 
The CK5Pro-GFP activity was quantitated using the sum of A) Fluorescence intensity and B) Area of GFP expression. The mean and SEM 
were obtained using four MCTS per treatment group and the data were normalized to the control group (ICQ nM P4 + vehicle). These data 
indicate that CSC phenotype is not downregulated as a result of cytotoxicity caused by 5FU and the estrogen receptor antagonist Tamoxifen. 



each individual plane as a distinct parameter. From the same 
image set, we can quantify the cross-sectional area of GFP 
positive cells, and the percent of spheroid cross-sectional 
area that express CK5Pro-GFP. For the compounds tested 
(RU-486, isotretinoin) we see concomitant decreases in 
overall CK5Pro-GFP expression and fraction of CKSPro- 
GFP positive cell area. 

In a high throughput screen where we can monitor these 
fields independently using a single image set, we may find 
compounds that can affect just one of these parameters. For 
example, one can imagine a compound that decreases the 
fraction of GFP positive cells (or GFP positive area) while 
maintaining the average expression levels in the GFP posi- 
tive area. In a secondary analysis of such a compound, a time 
course of expression could be done to see whether the effect 
is due to diffusion rate of a compound through the spheroid 
or to a particular effect on a biological pathway, although 
determination of the pathway could require a significant 



amount of downstream work. Regardless, the point of this 
theoretical example is that we can monitor many different 
parameters in live MCTS simultaneously, determine empiri- 
cally which parameters change in response to treatment with 
small molecules, and even look for differential effects of 
compounds on multiple parameters, using a single image set 
from each well of a multi well plate. Next, to further validate 
the CK5Pro-GFP activity as a faithful reporter of the CSC 
phenotype we treated MCTS with the clinically used breast 
cancer therapies, including: the cytotoxic agent five 
fluorouracil (5FU), and the estrogen receptor antagonist Ta- 
moxifen (Fig. 3) [40-42]. Both 5FU and Tamoxifen treat- 
ments did not induce down regulation of CK5Pro-GFP ex- 
pression, indicating that CK5Pro-GFP activity is a specific 
reporter for progesterone induced CSC phenotype expansion. 

Having determined satisfactory image acquisition meth- 
ods, we went on to write the analysis algorithm. The 
algorithm was written using a building-block approach in the 
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Fig. (4). A robust algorithm can capture botli local (single plane) data and global (whole spheroid) phenomena. Using a specific field 
within a well of a 96-well plate an algorithm for analysis of CK5Pro-GFP expression in MCTS can be used to capture data in a single plane or 
globally by using maximal projections. Panels 1-4 (bottom right) show the maximal projection of 8 Z-planes. Panel 1 shows the stained nuclei 
(Hoechst, blue) and the CK5Pro-GFP expression (green). Panel 2 shows the regions of the image tiled with spots (arbitrary colors). Panel 3 
shows spots selected by fluorescence intensity and size (selected spots in green, rejected spots in red). Panel 4 shows the entire spheroid re- 
gion across which high-content parameters were summed. 
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Harmony and Columbus software packages (Perkin Elmer). 
The first step in any image analysis algorithm is image seg- 
mentation, so that regions of interest can be divided, or seg- 
mented, according to various phenotypes. In a typical high 
content analysis algorithm, the initial image segmentation is 
done on stained nuclei, and current off-the-shelf analysis 
software will do a good job of segmenting individual cells 
based on the nuclear stain alone or nuclear stain coupled 
with a cytoplasmic stain. This proves to be problematic using 
a three-dimensional approach, because the cytoplasm of a 
cell is often in a different focal plane than the nucleus, lead- 
ing to very inconsistent and misleading image segmentation. 
We overcame this problem by taking a region-based ap- 
proach to analysis rather than an individual cell-based ap- 
proach (Fig. 4). Briefly, we use the analysis software to iden- 
tify regions of increased green fluorescence relative to back- 
ground, either in individual focal planes (Z-stacks), or using 
maximal-projection images, which combine the highest in- 
tensity values across all image planes. All green regions 
were tiled with spots of varying size, from individual pixels 
to cells or cell clusters with elevated CK5Pro-GFP expres- 
sion. Spots were filtered based on their fluorescence intensi- 
ties as well as their size. The CK5Pro-GFP expression thus 
selected is next summed across the spheroid area, and quan- 



tified both as fluorescence intensity values across the se- 
lected region(s), and as actual area of GFP expression (|xm^). 
This algorithm can easily be adapted for other fluorescent 
readouts, such as Red Fluorescent Protein and its various 
forms (tdTomato, mCherry), or other reporters besides CK5. 
Thus the algorithm, with minimal modifications, can be ap- 
plied to any number of MCTS model systems to monitor up- 
or down-regulation of phenotypic reporters in order to de- 
termine small-molecule efficacy in multi cellular structures. 
For example, the algorithm can be applied to show whether 
small molecules are having effects only on peripheral cells 
(outer cell layers) or on entire spheroids: In one modification 
to the algorithm (data not shown), we can separate the sphe- 
roid into concentric rings based on the nuclear stain of the 
entire spheroid, and monitor effects on GFP expression in 
concentric regions, from the center to the periphery. 

Finally, to demonstrate the robustness and utility of the 
MCTS system we have developed, we show that Z-prime 
values are greater than 0.5 for a number of output parameters 
from the analysis algorithm (Fig. 5). We also show that a 
dose-response curve can be generated using the same pa- 
rameters and a range of isotretinoin concentrations in live 
MCTS (Fig. 6). The IC50 values for isotretinoin (Fig. 6) vary 
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slightly depending on the parameter used (2-3 nM), but are 
consistent with our previously reported IC50 0.8 nM for 
blocking P4-mediated CK5Pro-GFP expression using the 
trans isomer (retinoic acid) in monolayer cell culture [38]. 

CONCLUSION 

The MCTS model represents a robust in vitro 3D tumor 
model that can be used to bridge the gap of information be- 
tween 2D-cell based assays and in vivo models of human 
tumors. Although the MCTS model is amenable to HTS drug 
discovery significant hurdles exist that have restricted its use 
to low-throughput cancer drug discovery. Nevertheless, 
MCTS HTS assays have been developed for endpoint detec- 
tion but yielding limited information. However, the MCTS 
model offers researchers a "systems biology" approach to 
cancer drug discovery miniaturized into a single assay en- 
compassing aspects of the tumor microenvironment and het- 
erotypic tumor cell culture. Despite the challenges that re- 
main with such systems, we demonstrate in this report a sig- 
nificant advance in achieving this using high-content image 
based algorithms to measure local and global complex phe- 
nomena in live MCTS. We show that acquiring data in this 
manner can be done in real time and is statistically robust for 
primary high-content screening cancer drug discovery and 
development. 
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